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i We present a calculation of the /3-term for D* D'j within a heavy-light chiral 

I quark model. Within the model, soft gluon effects in terms of the gluon condensate 

^ • with lowest dimension are included. Also, calculations of l/rric corrections are per- 

. formed. We find that the value of /3 is rather sensitive to the constituent quark mass 

in ! 

I compared to other quantities calculated within the same model. Also, to obtain a 

T-H ' value close to the experimental value, one has to choose a constituent light quark 

> : 

. mass larger than for other quantities studied in previous papers. For a light quark 

(N ■ 

mass in the range 250 to 300 MeV and a quark condensate in the range -(250 - 
' 270 MeV)^ we find the value (2.5 it 0.6) GeV^"*^. This value is in agreement with 

O ; 

■ 

iu: 

'. I. INTRODUCTION 

> ■ 

^ ■ Strong and electromagnetic interactions involving heavy and light mesons and photons 

may be described by heavy-light chiral Lagrangians Q]. Such Lagrangians are determined 
by chiral symmetry in the light sector in addition to heavy quark symmetry. However, the 
coupling constants of the terms in the Lagrangian are in general unknown. In this paper we 
focus on the so called /3-term for D* D^y, which is studied in . This term in the 

agrangian is also relevant for calculation of processes like —>■ 27 ^5j, and —>■ e^e~7 



the value of (3 extracted from experiment (2.7 it 0.2) GeV ^. 



Within chiral perturbation theory alone, the value of j3 is not determined. In this paper we 
will calculate the /3-term to order l/m^ within a Heavy-Light Chiral Quark Model (HLyQM) 



developed recently 



This model belongs to a class of quark loop models jsj, 0, 3] where 
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the quarks couple directly to the mesons at the scale of chiral symmetry breaking of order 
1 GeV. In contrast to most other versions of such models, our version also in corp orates soft 
gluon effects in terms of the gluon condensate with lowest dimension 

At quark level the Lagrangian of our model [7] includes the Lagrangian for Heavy Quark 
Effective Field Theory (HQEFT) jl?!. This means that l/nic corrections can be calculated 
within the model. It has been problematic topin down a precise value for /?, and chiral loop 
contributions to -D* ^ D'j are sizeable This may indicate that the chiral expansion 

will be problematic for /3. 

II. THE VALUE OF (3 EXTRACTED FROM EXPERIMENT 

The width of the D*^ meson has been measured very accurately : 

r(D*+) = (96 ± 26) keV (1) 
from this value of the width it is possible to extract a value for the 0*0-71 couplingjisl 

From (0) we find the following radiative decay rate : 

r(D*+ ^ = (1.54 ±0.56) keV , (3) 

where we have used the branching ratio Br(D*^ — > D'^'j) = (1.6 ± 0.4)%[l9|. To one-loop 
in the chiral expansion and to order 1 / rric we have [2, y] : 

r(z>--z>-7) = ^yM^[-k + |^ + f^^l\ (4) 

3Md* I 3 3 rUc 8n J 

where / is the bare pion decay constant / = 86 MeV[3|, and cu is the photon energy 
uj = (M|,, — Ml)) I {2M£)*). The term l/rric is due to photon emission from the c-quark. The 
bare coupling gj^ is related to Qa^^^- In the Heavy-Light Chiral Quark Model described 
to a certain extent in the next section, it is estimated to be 9^^-^ 0.57 ± 0.05. This estimate 
includes chiral corrections to one loop and I/ttlq corrections 7]. This value is very close to 
(0), in the following we will use the experimental value as the bare coupling. The quantity 
(3d is the effective, renormalized (5 including chiral corrections for the case when the D* 
and D mesons contain a d-quark j^: 

Pd = P \l + gA^ (^s^ + 3eK + ^er^A +3 (e^ + ^ (-^Sn + Sk + ^er^^)) \ , (5) 
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where we have used the notation : 

Eq. (0) includes vertex corrections for the /?-term and field renormalization. Combining 
equation (jS)) and (jH), we find two possible solutions for Pa '■ 

(3d = (0.7 ± 0.5) GeV^^ , /3d = (4.1 ± 0.2) GeV"^ , (7) 
where we have used mg = rric = 1.4 GeV. These values correspond to the bare values 

/3 = (0.7±0.4) GeV"^ , /3 = (2.7 ± 0.2) GeV'^ (8) 
Dropping the chiral corrections in (jH) (/3rf = /3) would give the result: 

/3 = -(0.1±0.3) GeV-^ , /3 = (3.0 ± 0.3) GeV^^ (9) 

The values in (jH)) and Q are consistent with those obtained in P]. Both the first values of 
P in equations (jH)) and Q are too low in order to agree with D*^ data 0. Thus, we are left 
with the second values of (jHl) and 



III. THE HEAVY LIGHT CHIRAL QUARK MODEL 

In this section we will give a short description of the Heavy-Light Chiral Quark Model 
(HLxQM) to be used in this paper. The model is based on Heavy Quark Effective Theory 
(HQEFT), which is a systematic expansion in l/mq (where tuq is rric in our case). The 
heavy quark field Q{x) is replaced with a "reduced" field, Qv{x), which is related to the full 
field the in following way: 

Q„(x) = P+e-™«'"^'Q(x) , (10) 

where P± are projection operators P± = (1 ± 7 • v)/2. The reduced field Qy annihilates 
heavy quarks. The Lagrangian for heavy quarks is : 

Wft = Q:tv-DQy + (^-Cm^ct ■ G + itD^)ls) Qv + 0( V) (n) 

where is the covariant derivative containing the gluon field (eventually also the photon 
field), and a ■ G = a^^G^j^^t"-, where a'^" = ■^[7'^, 7'^]/2, is the gluonic field tensors, and 
are the colour matrices. This chromo-magnetic term has a factor Gm which is one at tree 
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level, but slightly modified by perturbative QCD effects below mg. It has been calculated 
to next to leading order (NLO) 0,Q- Furthermore, {iD^)lg = CoiiDf - Cxiiv ■ Df. 
At tree level, Co = Ck = 1- Here, Cd is not modified bv perturbative QCD, while Ck is 
different from one due to perturbative QCD corrections j2l|. 
The Lagrangian for the HLxQM is 

■^HLxQM = -^^HQEFT + -^xQM + -^Int • (12) 

The first term is given in equation (fTT|) (Note however, that only soft gluons are considered 
to be included in (fT^ ) . The light quark sector is described by the Chiral Quark Model 
(xQM) j^, having a standard QCD term and a term describing interactions between quarks 
and (Goldstone) mesons: 

-^xQM = X b''{^D^. + V^, + 75^,.) -m]x- xM,X ■ (13) 

Here m is the (SU(3)-invariant) constituent light quark mass and x is the flavour rotated 
quark fields given by: Xl = C,^<1l and xr = ^iR^ where = {u,d,s) are the light quark 
fields. The left- and right-handed projections qi and qn transforms after SU{3)l and SU{3)r 
respectively. The quantity ^ is a 3 by 3 matrix containing the (would be) Goldstone octet 
(tt, K, rf). In terms of ^ the vector and axial vector fields and in (fT^ are given by: 

V,^\{ed,i + id,S}) , A,^-'-{i^d,i-id,e) , e = exp(^n//) (14) 

where / is the bare pion coupling, and H is a 3 by 3 matrix which contains the Goldstone 
bosons TT, K, 7] in the standard way. In ()13p the quantity Mq contains the current quark mass 
matrix M.q and and the Goldstone fields through ^: 

Mq = M^ + M^75 , where (15) 
M^^]^{eM\e+ iMqi) and M^^ -l(^t^t^t_^_^^^). (ig) 

The interaction between heavy meson fields and heavy quarks are described by the fol- 
lowing Lagrangian : 

^Int = —Gh [Xa^aQv + Qv HaXa] , (17) 

where Gh is a coupling constant, and Ha is the heavy meson field containing a spin zero 
and spin one boson (a is a SU(3) flavour index): 

Ha ^P+(P;7M-^n"75) , :^ = 7°(^a)V. (18) 



5 



The fields annihilates a heavy meson containing a heavy quark with velocity v. 

Integrating out the quarks by using by using J TT|) . |T3j) and (fT7)) . the lowest order chiral 
Lagrangian up to 0{mQ^) can be written as l2| 

£ = -Tr [TTatv ■ VbaHb] - qaTt {H^Hki^^^Al] , (19) 

where iV^a — '^^baD^ — and the axial coupling is of order 0.6. Eqs. (fT9|) and (fTl|l will 
be used for the chiral loop contributions. (Note that the eqs. (|T7jl and fTIH) both contain an 



additional term 



These are, however, irrelevant in the present paper). 
To obtain ()19|) from the HLxQM one encounters divergent loop integrals, which will be 
quadratic-, linear- and logarithmic divergent. Thus G'\j times linear combinations of these 
divergent integrals have to be put equal to 1 (the normalization) and respectively. For 
details, see the Appendix. Within our model, the values for the regularized versions of 
the quadratic-, linear-, and logarithmic divergent integrals are determined by the physical 
values of (gg), g^, and / respectively. The effective coupling Gh describing the interaction 
between the quarks and heavy mesons can be expressed in terms of m, /, g^-, and the mass 
splitting between the 1~ state and 0~ state. Using the equations ()A1|) . ()A3|) . ()A5|) . and ()A6|) 
in the Appendix, one finds the following relations between this mass-splitting and the gluon 
condensate via the chromomagnetic interaction in (lllj) [3] : 



W Ph _2m^ ^_ (7r + 2)^ 



where 



G') = ^^ , Gj, = ^p , ^^^^Cm(AJ , (20) 

n 717] p 71 



(1 + 3^^) + 3 
P - , pUH) = |mg(M,. - Mn). (21) 



24|. 



When calculating the soft gluon effects in terms of the gluon condensate, we follow 
The gluon condensate is obtained by the replacement : 

Soft gluons coupling to a heavy quark is suppressed by l/mg, since to leading order the 
vertex is proportional to v^v^G"''^" = 0, being the heavy quark velocity. We observe that 
the mass-splitting between H and H* sets the scale of the gluon condensate. Within the 
pure light sector, the gluon condensate contribution obtained by using has apriori an 
arbitrary value. One may assume that it has the same value as in QCD sum rules, or it has 
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to be fitted to some effect, in our case the mass splitting between pseudoscalar and vector 
heavy mesons. 

n 

The l/iTiQ corrections to the strong Lagrangian have been calculated in [7]. They may 
formally be put into spin dependent renormalization factors. This means that (jl9|l is still 
valid with the replacement H W = H {Zh)'^, where Zh and the renormalized (effective) 
coupling are defined as: 



Z^' = l + '-l^^^ , (23) 
9 A = ( 1 i^i- '^dj^£2) ) (5-1 - d^g2) , (24) 



where 



3 for 0^ I 1 for H*H coupling 

iM = S = < 

-1 for 1^ —1 for H*H* coupling 



and 



The chiral Lagrangian /5-term has the form 



(25) 



4m lo/T 16 m 

^ ^ ^ 12m 6 487r 16 m ' 

^ iCK-2n){^G')) , (27) 



64m^ 71 

9. = , = . (28) 

^ (vr + 2) 6m ' 2 ^ ^ 



IV. THE (3 TERM FOR D* 



Cp = ^Tr[HHa-FQ^], (29) 

where = {^"^Q^ + ^Q^'^)/2, Q = diag{— 2/3, —1/3, —1/3), and F is the electromagnetic 
filed tensor. The (3 term can be calculated in HLxQM, by considering the diagrams in 
figure n These diagrams will be calculated following 24], both for the gluon and the elec- 
tromagnetic field. The first diagram without gluons is logarithmically divergent. Performing 
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FIG. 1: Leading order diagrams contributing to /?. The double and single lines represent heavy 
and light quarks, respectively. The double dashed lines represent heavy mesons, and the wavy lines 
represent emission of soft gluons ending in vacuum to make gluon condensates. 

the calculation of the rest of the diagrams we used the algebraic program FORM We 
obtained the following expression : 



LO 



-AINJ2 + 



Nr. 



32 + 37r 



2 '^""-'^ ' Att 4m4 \^ 144 y ■ vr 
Using the relation ()A5|) for the logarithmic divergent integral I2, we obtain: 



LO 



1 + 



56 + 37r 



71 



(30) 



(31) 



2m2 [ Anp Apm? \ 144 
As seen from figurelH (3 depends strongly on the constituent light quark mass m, especially 
for m below 250 MeV. There is a partial cancellation between large terms in ()3H1 . and for 
values of m smaller than about 200 MeV, (3lo turns negative. Note also that if the gluon 
condensate is dropped in , [3 would be too big, as seen from figure El 

Apriori one might hope that l/rriQ corrections stabilizes (3 for values of m in the range 
190 to 250 MeV used in 3]. However, as seen in the following, this will not be the case. The 
l/mg corrections to /? can be found by calculating the diagrams in figures |21 and El We 
find that the kinetic term in (jllj) give a negative contribution to (3, and the chromomagnetic 
term a positive contribution. However, the kinetic term dominates in absolute value, and in 
total, l/niQ corrections give a negative contribution to /3 for m below 250 MeV. 




FIG. 3: l/mg diagrams from the chromo magnetic operator, contributing to /? 
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FIG. 4: (3lo as a function of the constituent light quark mass for parameters from ref. 
Combining the leading order and the 1 / uiq corrections we find : 



(3- 



ZhZe'I^lo + 



where 



-4(1-^^)- 
(7r + 4) 



4m 



2f - 



27r2 



+ 



1 

12 



-C 



A/- 



72 



,Ol.s 



m? vr 



G') -CK[f 



TT 



192 



,a 



(32) 



where rg contains the result of the diagrams in the figures El and El A plot of (3 is shown 
in figure jSl As we see, the result varies significantly with m. This is in contrast to other 



quantities studied in 
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ll|. We observe that the quantities studied in |7,lll| have contribu- 



tions of zeroth order in m. In contrast, (3 is of order 1/m, which means that (3 is expected 
to be more sensitive to variations of m. The l/mq corrections are not sizeable, even if the 
charm quark is not very heavy. However, they pull in the wrong direction compared to the 
experimental value of order 2-3 GeV^^. To obtain a value closer to the experimental value 
for /5, we need a value for m higher than used in 0, However, this will lead to values 
of fs and fr, which are too small . This can be compensated by using a higher value of 
the quark condensate iqq) than the one obtained in QCD sum rules, which was the one 
used in p, |11|| . This might be acceptable because it is not clear that our model dependent 
quantities {qq) and (^G^) should be exactly those obtained in QCD sum rules. Taking 
a higher quark condensate and making a new fit we obtain the result given in the tables I 
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FIG. 5: /? as a function of the constituent light quark mass for gji^ = 0.59 

and II. Note that the curve for (5 itself as a function of m does not change, but we are now 
allowed to go to higher values for m at the curve. 

Chiral Lagrangian terms with extra current quark mass can be obtained by taking the 
derivative of the expression for (3 in (jHUj) with respect to m, when keeping Gh and (^G^) 
fixed. This gives: 



ea 



where 



FM 



a 



Tr 



Gff 



Nr /32 + 37r\ ,cts^2\ 



+ 



Am I vr^ ' V 72m^ / tt 



(33) 



(34) 



This term is unimportant for D* D'-f because of the small current quark masses for the 
u, d quarks, but it gives a sizeable contribution to -Ds7. 

V. CONCLUSION 



We have calculated the quantity (3 and found that it depends significantly on the con- 
stituent light quark mass m. To obtain a value close to the experimental one has to pick 



111, where 



values on the high side (m= 250 to 300 MeV, say) compared to the values in l7, 

un 

m = 220 ± 30 MeV were used. We have redone the fits in [21 at the price of a higher 
value of the quark condensate {qq) in the range —(250 to 270 MeV)''. The results are given 
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TABLE I: Predictions of HL^QM for different input parameters in the L'-sector. 



{qq) = -(230 - 250)^ MeV 
m = (190 - 250) MeV 


(qq) = -(250 - 270)3 MeV 
m = (250 - 300) MeV 




I / zn u.u J vjre v 




(6.4 ±0.4) GeV-i/2 




I'oqn _|_ on\ MpV 




(330 ± 20) MeV 


9i 


I'n o _|_ n 1 \ rjpV 
I u.o in u.± ) vjc V 




(0.61 ±0.07) GeV 


y2 


I'n Qo _|_ n n4.\ rjpV 




(0.25 ±0.02) GeV 


£i 


_/'n o _L n q\ p^v 

— l^U.O ZC kJ.OJ vjt; V 




-(0.6 ±0.2) GeV 


x. 


8-1-09 
U.O It U.Z 




0.5 ±0.1 




(0.34 ± 0.05) GeV^ 




(0.29 ±0.04) GeV^ 




(220 ±60) MeV 




(220 ±35) MeV 


fD* 


(260 ±85) MeV 




(235 ± 50) MeV 




(245 ±70) MeV 




(260 ±45) MeV 


Id* 


(280 ±95) MeV 




(270 ±60) MeV 


Id* I fo 


1.2 ±0.1 




1.12 ±0.04 


fos 1 Id 


1.18 ±0.06 




1.26 ±0.04 


P 


(0 ± 2) GeV^i 




(2.5 ± 0.6) GeV-i 



in the tables I and II. For the case D* —>■ D'-f we have found that (the bare) (3 = (2.5 ± 0.6) 
GeV~^ to be compared with (3 = (2.7 ± 0.20) GeV~^ extracted from experiment. 
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APPENDIX A: THE LIMIT m ^ 

In this Appendix we will discuss the limit of restauration of chiral symmetry, i.e. the 
limit m 0. In order to do this, we have to consider the various constraints obtained when 
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TABLE II: Predictions of HL^QM for different input parameters in the i?-sector. 







g^D*Dn ^ Q_5g ^ 0.08 fl8j 




{qq) = -(230 - 250)3 MeV 


{qq) = -(250 - 270)3 MeV 




m = (190 - 250) MeV 


m = (250 - 300) MeV 


Gb 


(8.3 ±0.7) GeV-i/2 


(7.2 ±0.5) GeV-i/2 


\ -n ^ 1 


(300 ± 25) MeV 


(340 ± 20) MeV 


9i 


(1.4 ±0.3) GeV 


(1.0 ±0.2) GeV 


92 


(0.39 ±0.05) GeV 


(0.31 ±0.03) GeV 


El 


-(0.9 ±0.4) GeV 


-(0.6 ±0.2) GeV 


Ai 


1.1 ±0.3 


0.8 ±0.1 




(0.42 ±0.03) GeV2 


(0.39 ±0.03) GeV^ 


Jb 


(190 ±50) MeV 


(185 ±30) MeV 




(200 ±60) MeV 


(190 ±35) MeV 


fB 


(210 ± 70) MeV 


(215 ± 45) MeV 


fB' 


(220 ± 70) MeV 


(215 ± 45) MeV 


Ib* I Ib 


1.07 ±0.04 


1.02 ±0.02 


IBs I Ib 


1.14 ±0.07 


1.22 ±0.02 


^Bd 


1.51 ±0.09 


1.52 ±0.07 


Bb. 


1.4 ±0.1 


1.4 ±0.1 




1.08 ±0.07 


1.16 ±0.04 




-(2 ±3) GeV-i 


(1.2 ± 0.8) GeV-i 



constructing the HL^QM [3]. 

To obtain (fTT^ from the HLxQM one encounters divergent loop integrals, which will in 
general be quadratic-, linear- and logarithmic divergent. For the kinetic term we obtain the 
identification: 

-^GlN, (ly, ± 2ml, + ^^^{^G^)] = 1 , (Al) 



and for the axial vector term 



-iG%Nc — h/2 H h 2m/2 ± I- ^ — G^) \ = qa ■ (A2) 
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Combining ()A1|1 and ()A2|) . the strong axial coupling g^^ niay be written: 

4 / ZTTl 

9a = ^ ~ ^Qa , where Sqa = - ^iOj^N^ i I3/2 - ) • (A3) 



Here the term 1 for (7^ corresponds to the normalization ()A1|) and the term 6q A is a dynamical 



deviation from this. It should be noted that also within other models j22, |23[ , may be 
written as (7^ = 1 — 6gj(, but the expression for 6gj( is model dependent . We observe 
that the (formally) linear divergent integral I3/2 is related to the strong axial coupling 
Analogously, within the pure light quark sector (the xQM), it is well known that the 
quadratic and logarithmic divergent integrals are related to the quark condensate and /, 
respectively ^, 
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14 



3 



(qq) = -AtmNJ, - -^(^G^) , (A4) 

12m TT 

P =-iAm'NJ, + -^^(^G') . (A5) 

Eliminating 13/2 from the eqs. (jAlll and ()A3|1 and inserting the expression for I2 obtained 
from ()A5|) we find the following expression for Gh '■ 



G%, = '^dA) where 771 = — . (A6) 

2 f ^ + "^^^ — (^Q^\ 32 



The divergent integrals Ji, I2 and 13/2 are (for = 1, 2) 



/ d'^k 1 ^ f ji^ 1 



where J3/2 is proportional to the cut-off in primitive cut-off regularization: 

^3/.-4(l + 0(^)) . (A8) 

where the cut-off A (which we only have used in the qualitative considerations in this Ap- 
pendix) is of the same order as the chiral symmetry breaking scale A^. In contrast, 13/2 is 
finite and proportional to m in dimensional regularization. Within our model, the values 
for the regularized integrals Ji, I2 and 13/2 are determined by the numerical values of ( gg ), 
/ and (yfyi, respectively. 

Looking at the equations ()A4jl and ()A5j) . one may worry j26^ that {'qq) and / behaves 
like 1/m in the limit m — unless one assumes that (^G^) also go to zero in this limit. We 
should stress that the exact hmit m = cannot be taken because our loop integrals will then 
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be meaningless. Still we may let m approach zero without going to this exact limit. In the 
pure light sector (at least when vector mesons are not included) there are no restrictions on 
how (^G^) might go to zero. In the heavy light sector we have in addition to ()A4|) and ()A5|) 
also the relations (jAlj) . ()A3|1 . and ()A6j) which put restrictions on the behavior of the gluon 
condensate {^G^) for small masses. As (^G^) has dimension mass to the fourth power, 
we find that (qq) and may go to zero if (^G^) goes to zero as m'^ or m^A (eventually 
combined with ln(m/A)). However, the behavior m^A is inconsistent with the additional 
equations ^ and Still, from all equations (jUI) and (jHl), we find the 

possible solution 

a 1 

{—G^) = cN^m^K{m) , where K{ni) = {-Aih + — ) , (A9) 

and c is some constant. Then we must have the following behavior for G^, and when 
m approaches zero: 

1 77? 7T7^ 

Gh-j^ , (1 + 3^^) ~ -K{m) , ^^l ~ —K{m) , (AlO) 

with some restrictions on the proportionality factors. Here, the regularized I2 is such that 
for small m, K{m) = (ci + C2lnm/A), ci and cobeing constants. The behavior of G^ is 
in agreement with Nambu-Jona-Lasinio models 9|. Note that in our model, 5g_/( 4/3 



22 1 for a free Dirac 



(corresponding to — — 1/3) for m — 0, in contrast to dg^i 2/3 in 
particle with m = 0. Note that in [7] we gave the variation of the gluon condensate with m 
for a fixed value of fiQ. For the considerations in this Appendix, we have to let fiQ go to zero 
with m in order to be consistent. When m — 0, we also find that (3 1/A, provided that 
the coefficient c in ()A9|1 is fixed to a specific value ( which is c = 576/(37r + 32) (1.93)'*). 

Before closing this section, we will shortly comment that in the limit where only the 
leading logarithmic integral I2 is kept [3], when gA 1 lli] and p — > 1, we obtain from 



equation (jHT|) the non relativistic quark model result 



^ /3nr = - . (All) 
m 
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